
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



ON THE NABLA OF QUATERNIONS. 

By Mb. Shunkiohi Kimuba, Tokio, Japan. 

In studying McAulay's differential linear vector function, I came across 
the nablas whose arguments are the sums of vectors, which led me to consider 
the same operator in its wider sense, thereby extending it to the case of 
quaternion argument with transformations that might be useful in attacking 
the problems of physics and geometry. The discussion is not exhaustive, but 
it is hoped that the prominent features of the operator have been touched 
upon. Along with the quaternion treatment are appended the results in semi- 
Cartesian forms, thereby showing how much simplicity is gained by the former 
treatment. 

Throughout the paper, p, a, q, p have the following meanings whenever 
they are used in semi-Cartesian treatment : — 

p = ix+jy + ks, 

<r = ix + jr+ kz, 

q = t + ix+jy + kz, 

P = T + ix+jr + kz. 

Occasionally the same symbols are used with different meanings in different 
divisions of the paper ; these will, however, be carefully denned or explained 
in each case. 

§1- 

In stating the differential symbols in the original definition of a nabla, 
two kinds of notations are used ; thus, for example, 

^ = l 9i +J 9y + k 9z'\ 

> (1) 

. d , . d . ; d | 

These are immaterial when p is the independent variable vector, i. e. when *, 
y, z are the three independent scalar variables ; but when we come to gen- 
eralize the same operator we are restricted to one of these two notations. 
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Thus, for reasons to be seen later on, I propose to write down for any vectors 
p and a, 

. 9 , . 9 , . 9 ■) 



9x 

3 



3y 



I 



*V ==z dX + J^' + k 



3 
"9Y 



9Z' ) 



(A) 



30 

where -^ , etc., represent partial differential coefficients of a certain function Q. 
Quite reasonably we may extend the nabla in the following form : — 



9 , • 9 , . 9 , 7 9 



kqV 



9_ 
3t 



.3 .9 
J Sy 



(B) 



9x 



'3z- J 



When q is a function of t, or when it is an independent variable quaternion, 
that is, when t, x, y, z are the four independent variable scalars, then invari- 
ably we have 

d — SdKq a v = Sdq A - a v 5 (2) 

and this is the form usually adopted for the case when q degenerates into a 
vector. 

If, on the other hand, we were to express a ^, g y in the form of total differ- 
ential coefficients, 

d , ■ d , . d , 7 . d 
~di 



«V = ^h + i ^,+Ji;, + k jz> 



dx 



dy 



dz> 



we could obtain the form (2) only when t, x, y, z are independent ; if q were a 
function of t, we should have 

SdK ** = dt \di + [di\ 9x+[l\ 9y + [dij 9z\' 

while by definitions (A) and (B), even when q is a function of t, 
d = SdKq q v 



f 9 .dx 9 .dy 9 dz 9 1 
\Wt ~ t ~ di 3x + di dy + dt 9z j 



(3) 



As illustrations of the advantage of using definitions (A) and (B) may be 
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cited the operator employed in hydrodynamics, which is obtained in the form 

d_9dx9dy9dz9 
di ~ 9i di 9y + ~di 9y + Tt 9s 

and Euler's equations, which are written at once in the form, 

( Vdq\ S d fi aVl = VqV Qdt . (5)* 

Thus it is seen that nablas in their extended form have direct physical appli- 
cations, not to mention the " electrodynamic quaternion " of Hamilton. 
When q degenerates into a scalar t, to be consistent we must write 

X = j t , (6) 

and when t is the only independent variable, it is immaterial whether this be 
written in the form of a partial or a total differential symbol. 

§2. 

Usually the differential operator is written in the following forms, when 
q, p, a, and t occurring therein are each taken for the sole independent vari- 
able ; and these are taken as new definitions of nablas : — 

d = SdKq q y = Sdq K( tf 

= SdKptf = Sdpxrf 

(7) 
= SdlTa^ = Sda Ka ^ 

= SdKtfi = Sdt K fi , 

where it is to be noticed that when q degenerates into a vector or a scalar, 

Kp = — p , ^v = — „V = K ?V ; ' 

Ka = —a, j^v = — „V — ^<rV 5 \ (8) 

Kt = t , ^-(V = <v ; 

and, in general, 

KqV = KqV • 

* Kirchhoff, Mechanik, S. 164. The subscripts, 1, show that Vdg is the quantity to which q /\ 
applies. 
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Similarly with respect to Laplace's operator, we usually call y 2 the nega- 
tive of the operator, but since 



c o & , 3* P 2 , 9 2 1 



5* 2 



Sar 1 ' 9tf n 3s 2 



we have, as may readily be seen, 

8qVqV = SxqVztV = 



2 2 



3i! 2 3* 2 
In reality, the positive of the Laplace's operator is 

§3. 






3s 2 



(9) 



(10) 



For the change of the independent vector or quaternion variables, if we 

take a to be a vector function of p, and p to be a quaternion function of q, we 

have always 

da = tpdp , dp = <Pdq , } 

> (11) 

dp = <p~ l da ; dq = 0~ J dp , } 

where <p, <P are respectively a linear vector and a linear quaternion function. 
Then, since 

d — Sdp K ^ — Sf-^dffxtf = Sdof'- 1 ^ 

= Sdo Ka y = S<pdp Ka A = Sdp<p' K ^ , 
and 

d = Sdq Kq y = S4>- x dp Kq v = SdpV- 1 ^ 

= Sdp x rf = S<Pdq Kl & = Sdq<D' Kl x , 
we obtain the relations 



pV = <P cV , 
k p V = ®'~ 1 k<,V = <F~ l K q V or 



<N = <P'pV > 
P V= «>'-\v 



(12)* 



(13) 



It is almost needless to remark that the nablas of quaternion, vector, and 

* Tait, Quaternions, 3rd ed., p. 373 ; <p', <t>' being conjugates of <p, 4>. 
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scalar arguments are respectively quaternion, vector, and scalar operators, of 
which the first is the sum of the last two : — 



V«V = — Vk q V = — V ^qV = v„V = pV • J 
The transformation of Laplace's operator is effected by 



(14) 



> 



and 



(15) 



where <p<p' and <?^ are self conjugate functions. 
Thus, for example, when 

a = Wp , 

where Wis a scalar function of p, then 

da = — W l pS^idp + Wdp , 

in which W r is not a factor, but the quantity to which the operation y, applies, 
which is shown by the common suffix.* Hence 

<p X = WX — W lP S Vl X , 

<p'X = WX — v WSpX , 

0X = W 2 X - WV( V WVpX) , 

Stfhpp : SXp = W\ W — Sp V W) , 

where X is any quaternion, scalar, or vector, <p is Hamilton's inverse function, 
and v is used for p y when it is already operated, and we have 

w„V = <P'~\V = ( W fV + V WSptf - pV Sp V W) :W(W- Sp V W), ) 
P V = f'w P V — WVpV — V WSp Wl tf , V (16) 

^VrpV = - "^>rpV>rpV + ™S V W w , v Sp Wl>V ~ (v WVS*p w # . J 

The case a = p n (n = 2m -j- 1) is included in the above, for here we have 

W= (— l) m (TpY m , 
* See McAulay, Utility of Quaternions in Physics, p. 13. 
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where Wis a constant, equal to n ; whence y TP"=: 0, and 

1 

npV = - pV • 

§4. 

The results of the above transformations may be expressed in a different 
form without the use of <p ; for we may write, in general, 

da = — i (afi^dp + a 2 S\?. 2 dp + a 3 ^ 3 dp) . 

Hence f of the last section is of the form, 

<F = — 4 K#Vi r + *2#W + *Av) • (17) 

Whence, the third degree invariant m and the conjugate of Hamilton's inverse 
function of <p are 

m = \ 8a x a 2 a % ^ViVzVs . 

fir = 6 ( y°\ a % ^vWi*- + Va 2 a B ^VsW + Va i a l S^^t) 

= i V^i<r 2 ^v 2 v^ ; 

also 

tp'T = —\ (Vitfffir + y t Sa 2 T + Vs&V) 

= — ^ l Sa,T . 
Hence we obtain the following results of transformation : — 

P V = <p'*V = — Vi ^wV > 
ffV = jp'-Vv = m-^'pV 

= 3 [ F^ /SpWiVzl '• C&Ws <Sv'iV»V»] . 
SpVzeV — — ^ViVz S ff wV 8 a wV • 



(18) 



In these expressions the suffixes applied to a and v are due to McAulay 
and mean that the differential operator v applies to a of the same suffix while 
as a vector y remains in its place, y is used for p y when the same is already 
applied to vectors a in the expressions, while P v remains in its nascent state. 
The brackets in the expression for ff v indicate that the dividend and divisor 
are independent with respect to the suffixes, 1, 2, occurring in both. 
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The semi-Cartesian forms for the above expressions are as follows : — 
Vofa /S'v2Vi r = — 3 </>'t 

„ f Tr So So m , T7 -So So a- i rrSo So c,. 



3x Sy 



3o 



% ySz 



Sz dx 



3a 



o lSVl r = _ Sfr = £ «r + = Sjr + ^Skr, 



ViSo^t — — B<p'z 



Sx Sy Sz 



3s : 



(19) 



and so on. Transformations of quarternion nablas are treated in the same 
way. 

§5. 

If we try to deduce the results of the last paragraph from purely Car- 
tesian forms, we arrive at somewhat novel but complicated expressions. Thus, 
supposing o to be a function of p, and consequently vice versa, we have as 
usual 

dX = - 7 — dx -f ^— dy + -^- dz , 
9x ^ Sy * Sz ' 



iXr SY , , SY ■, , 3Y , 
dY = -=— dx + -tr—dy + ~x—dz, \ 

8x 3y " z>* < 



dz 



(20; 



3x 3y y Sz 

SolviDg this system of equations we obtain 

M a dx = M Xx dX + M yx dY+ M Zx dZ , "j 

M„dy = M Xy dX + M Yy dY+ M^dZ, \ 

M a dz = M Xz dX + M Yz d Y + M Zl dZ , J 

where M a is the determinant 

9X 9X 3X 
~Sx ' ~3y ' ~~Sz 



(21) 



M. 



SY 3Y_ SY 
3x ' 3y ' Sz 

SZ SZ SZ 



(22) 
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and Mx„ etc., are its minors. Consequently 

3x _M Xx 3x __M Yx 3x 
3X~ M„' 3Y~~W a ' 32'' 






3X~ M> 3Y~ M' 3Z~TI' 



3z 
3X' 



M Y 



3z 
3Y 



M a 



3z_ 
3Z : 



M,, 



M„ ' 3Y~ M a ' 3Z~ M L 
In quarternion notations these take the forms, 

M a = ■kSa l <r 2 ff ;i S'v 1 'v 2 'v i , 

M 3Y3Z 3Y3Z , „. c . 



M, 



Xy 



3Y3Z 
3z 9x 



9Y3Z 



9x 8z 



= \Sia x a 2 aS^Vs 



M 3Y3Z 3Y3Z , Q . Q , 
M ** = -3x-3v;-^y~3-x- == iSl<T ^ ^ VlV * 

irf 3Z3X 3Z3X . . . 
My * = 3J -W ~ & ^T = ^^ &VlV2 



^ PZS^T 3Z3X . . „. 

Mv " = -3z-~W~3x~-3z- = iSja ^ ^ VlV * ' 



L Yy 



M _3Z3X _3Z3X 
Y " 3x 3y 3y 3x 



= h^J a \ a 2 8&ViVt > 



m 3X3Y 3X3Y . OJ a . 
Mz * = Jy--3z-~^'3y- = * 8ka ^ Sl ™* , 

w 3X3Y 3X3Y xVU „. 
Mz » == ~3z-'3x--^-3z == ****>** S - ? ™* ' 

3X3Y 3X3Y 



M„ = ^^L 



3x 3y 3y 3x 



= ^Ska^o* Sky{y 2 



(23) 



(24) 



Thus we may observe that there are at most nine independent quantities, hence 
these can be taken together in a linear vector function. By employing these 
quaternion forms in the expression, 

3 , . 3 . , 3 



V = * 



3X 



+ JjY + k 



3Z 



si [2e 2. jl 2m A 4. A 2.1 

\_3X 3x "•" 3X 3y ^ 3X 3z\' 
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we obtain the same results as in (18). The same kind of treatment applies to 
quaternion transformations. 

§6. 

Expressions of partial differential coefficients in terms of ip of § 3 lead 
to many interesting results. 

Since by (12) and (13) we obtain 



^y = — Si„v — — xSpVf l i , 



m7 = - Sj,v = — #pW M > f 



Yi = — *^w = — ^.Vf lJc ; 



^ = — Si P v — — S^fi , 
= — SjpV = — S °V<PJ , r 



ex 

d_ 
3y 

3_ 

3z 



Sk a 



S a Tj<pk ; 



3 



fr = - s Jpv = - s^<p- i j, 



3 
ft 



3x' 

3_ 
3y 

9 

3z 



— Sirf = — S p y0i , 

— Sj„V = — S pV@J > 

— Sk^ = — S p Tj0k ; 



(25) 



(26) 



(27) 



(28) 
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we obtain for the transformation of vector variables, 



= — Sif-H , JL=— Sjf-H , jJ^ = — Sk<p-H ; 



sX 

Sx 



S~X 

Sz 



jr = -Sif- l j, ^L = -Sj f - l j, Jr 



S& f - l j; 



^=-Si^k, %> z = -Sj<p-k, ^ = -Sk^k; 



> (29) 



9 -X=-Sifi, d -*- = -Sj<fi, 3 ^=-Sk f i; 

SX . . SY a . . SZ cij • 

-fy=- St W> -ty^- S M> fy^- Sk W> 

9 -§ = -Si<fk, 9 -^=-8j f k, 8 -Z=-8hpk; 



(30) 



and for the transformation of quaternion variables 

i<_ = SV~ l i. -^L = S<P'-H H-= so?- 1 ! — = 8#- l k- 



it = - ^-h- , J5- = - #**-% # = - ^'^ , 41 - - -at*-^ ; 



5JT 



£T 



3X 



5X" 



^ = - S0~\j , l x T = - SW-y, ^=- Sj4>->j , *, = - ^(P-V; 



> (31) 



§ = - 8+*k, -* = - SW-% H = - ##->*, £ = _ aw*. 



St ' St ' St J ' St ' 



ST 

3x 



... sx 
' Sx 



SY 
Sx 



%± = — S0i , ^ = — SiM , ^L = — Sj0i, ~ = - Sk0i ; 



sz 

Sx 



' ;r = -S<Pk,~ = - Si<Pk, -~=- Sj0k, ~ = - Sk<Pk. 

<vV. W«*. £*JJ 



S (32) 



3s 



Pz 
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Thus, when <p and are self-conjugate functions, we obtain 

V.vp = . ) 
V pV a = 0-S 

V„VP = o , 

r pVq = o. 

Expressed in Cartesian symbols these become 



3Y _3X 
3x 3y ' 


3Z _ 
3y 


3Y 
' 3z ' 


9X_3Z 
dz 3x 


3y 3x 

3X ~ IT' 


3z 
3Y~ 


3Z' 


9x 3z 
3Z~3X 
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(33) 



(34) 



We shall notice here that the second equation of (33) and the first of (34) 
are respectively a condition of integrability of 



Sodp = ,* 
Spdq = , 



(35) 



in the forms fp = c, Fq = c, respectively ; and (33) and (34) show that when 
these are integrable then Spda = and Sqdp = are also integrable, which 
is also at once seen from the identities 

dSpo = Spda -f Sadp , 

dSpq = Spdq -\- Sqdp . 

Again, if <p, are self -conjugate, we have for any constant vector a and 
any constant quaternion a, 



(36) 



^Sip = — <p H , 


tfSio = — (fi , 


*v s Jp = — <r l J . 


pV Sjff = — <pj, 


a ySkp = — (p~ l k ; 


p ySk<j = — <pk ; 


a ySap ■= — <p~ l a : 


rfSao = — <pa : 



* Tait's Quaternions, 3d ed., p. 253 ; Hamilton's Elements, Art. 415. 
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p ySq = K-qSQ-% gV Sp = KyjSQrj , 

pV Sig = — 0-H , qV Sip = — 0i , 

pvSjq = - *-]j , q vSjp = -&j ,y 



(37) 



p ^Skq = — <P- X k ; tfSkp = — <2>k ; 
p TjSaq — 0~ l Ka : q ySap = — <PKa: _ 
in which the -y-pairs are extensions of £-pairs, such that 

Q 0?, i) = Q («vi» ft) = # (ft, «vj 1 

= Q (1, 1) + Q (i, t) + # (i,i) + £ (*, *) I (38) 

= <? (1, 1) + Q (C, C) , J 

£ (ify 7J)=Q ft, Krj) = Q Uv.> ft) 

= # (ft, K Q V) = Q (Kq» , V ) = £ ( aV , ^) J. (39) 

= 0(i,i)- «(C,0; 

where $ is any quaternion function linear in each of the arguments mentioned. 
Now, whatever the linear vector function <p may be, we have 

3a 



dp 
JX 

dT 

dp 

'S'Z 



<p 



? y 



dx 



(fi 



da 

* = " 



3 P= f -'k. 



3a 

3z 



= <pk. 



(40) 



Hence, we obtain 






TO 
TO 



rt a = C^C = 2 r — ™" , 



<5pV ff : 



TO 



T 7 pV ^=2 r ; 



(41) 



(42) 



where m, to', to." are the three invariants of <p, viz.: 

<p s — m"<p 2 -\- m'<p — to = 
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and y, 3 are rotation vectors of <p and its Hamiltonian inverse function <j> ; thus 
m = — lS<p'ZK> m' = — 8&Z , m-=—SCfCO 

Similarly whatever the linear quaternion function may be, we have 



(43) 



31 



dq 



St 



ST J ' Sy 



K9i, 



K0j, 



(45) 



Thus we obtain 



&—****■■ £— **j 






(46) 



(47) 



Thus by (34), (46), (47) we have for any self-conjugate linear quaternion func- 
tion 0, 

yr V0* = , Vc V0~X = . (48) 

This means also, that when is self-conjugate, the linear vector function cp 
derived* from it is also self-conjugate, where <p is defined for any vector r and 
any quaternion s, by 

if Vs = V0 Vs ) 
or > (49) 

ipr = V0r . ) 

Hence in (31) and (32), if «, ft stand for either of the unit vectors i, j, k we 
have 

Sa0p = Satpp , Sa0-'p = Sa<p- l p . (50) 

* In regard to the reduction of <p into <p, see Hamilton's Elements, Art. 364. 



140 KIMOBA. ON THE NABLA OF QUATERNIONS. 

§7. 
In the definitions of p, q, p y, „y, if we put 

p x = ix, p 2 = jy , p s =• kz ; 
q l = t, q 2 = ix, q 3 = jy , q t = kz ; 

and further take t, x, y, z to be four independent variables, then the p's, and 
^'s are entirely independent of each other both in tensors and in versors ; and 
we have 

P V = 2 Pl V = ^'p,V » ffV = s«,V = AiV • 

Consistent with this, if p, a, r, . . . are independent vector variables, and 
if p, q, r, . . . are independent quaternion variables, we may write 

(51) 







2pV 


' = ^pV : 


> 


IqV 


= ^«v- 


However, 


although 


we may 


write 












p,V 


= *»V = 


. 9 
dx 


= 


TT d 

Up "JTp- x > 


we cannot extend this into 














pV = 


TT 3 

Up JTp> 


or 


«v 


TT d 

= Uq JTq> 



since in the first case Up x is a constant vector, while in the second Up or Uq 
is a variable vector or quaternion. Now to prove (51) we have merely to 
write 

d = SdKIp SpV = ISdKp pV = SdKIpS^ , 
and 

d = SdKSq XqV = ISdKq qV = SdKSqI qV , 

in which, on account of the independency of the constituents, in differentiating 
any function Q, such terms as 

SdKp^ or SdKq p -y 
do not exist. 

However, when the p's or ^'s are functions of one of themselves, or when 
they are functions of some other vector or quaternion, such simple results as 
(51) cannot be obtained ; for, since in this case 

d2p = (2y) dp , 

dZq = (10) dq , 
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we have 



(52) 



seV = (-2?)' VV . or pV = (^V)' SpV ! ) 

s,V = (^tfy-'.V . or , V = W*»V- > 
By the definitions (A), (B), (C), whether t, x, y, z are independent or not, 



we have 



KtV f = tVKt = 


i, 


K„W = pV K f> = 


3, 


Km = «v K Q = 


4, 


m = jkV-Ks = - 


-2. 



(53) 



But, if we take n quantities, />, a, r, . . . , and n other quantities, jp, #, r. . 
to be independent vectors and quaternions, we have 

z t \r2't = 21 = n 



■ztflKp = 


2'd = 


xqvZKq = 


14: = 


■sJtSq = 


-12 = 



3n, 
4n , 
2w. 



> 



(54) 



write 

we have also 



y 



(55) 



Since, however, whether the constituents are independent or not, we may 
lt = T, I P = P, Iq = Q; 

%ty2t = T ^T = 1, 
ifi/SKp = rV KP = 3, 
i^IKq = (fiKQ = 4, 

jtfl'q = Q yQ = — 2. 

Hence we must distinguish between the two cases, in the first of which the 
constituents of the arguments are taken separately, while in the latter they are 
taken conjointly, and for the latter case I propose the notations, 

(2t)V » (Sp)V » (2«)V • (56) 

In connection with this, I might mention that when we take nine scalars 
of any linear vector function <p as independent variables, then, as McAulay* 
has shown, 

iQ=-QiSdfC 4 g l c. 

•Utility of Quaternions, etc., Art. 5. 
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This may also be written 

That is, the total variation of Q is equal to the sum of the variations of the 
same due to the variations of three vectors (pi, <pj, <pk which are independent of 
each other. 

§8. 

Seeing that the nabla of a quaternion argument is capable of direct phy- 
sical application, I propose now to transform such a nabla into the same 
operator with a vector for its argument, and vice versa. 

The simplest case is when the vector element of the quaternion argument 
is the vector under consideration, thus 

«V = sqV + vqV = s«V + pV • (57) 

In general, let p be a quaternion function of a vector p, and let a be a vector 
function of a quaternion q ; then we have 

dp = 4>dp , ) 
da = VWdq ; 

where 0, W are linear quaternion functions, having no relations between them. 
We cannot produce the second from the first, nor the first from the second ; 
for 01 is not given in the first case, nor SWdq in the second. 
Thus, 

d = Sdp Kp y = S0dp Xp y = Sdp0 r Kp y — Sdp Xp v 

d = Sdo Ka y = S¥dq Ka .y = Sdq W' Ka y = 8dq Kg y ; 



I (58) 

; > 



whence 



or, also, 



^V = V<P KpV _, or , v =:-rVK P vi 






(59) 



pV = fVi'8Pi x P V > «V = «Vi 8 ff i k«V ■ (59') 

Thus Laplace's operators become 

StfKtf = S VVK pV VVK pV = ~S V0- KpV VV KpV , 1 
S Q y KQ y^SW' aV KW aV = S <rV WKW' a y, v (60) 



(61) 
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Since the general form of a linear quaternion function is 

4>q = la'Saq — a' Sag + b'Sbq -+- dScq + d'Sdq , 
we have 

<% = laSa'q , K4>q = IKa'Saq , 

KWq = ZKaSa'q , 4>Kq = Za'SaKq , 

(K$)'q = laSKa'q , K$q = IKa'Saq , 
and so on. 

Thus when q, a are functions of t only, we write t \r for -^- since it is the 

sole independent scalar variable ; then 

*V = «Vi ^ v i k„V = tVi 8 w 1\ k*V = tVi Sq { Wk,V = tVi $2i jt« V- 
Again, whenj? and p are functions of t only, 

These results agree with those obtained before. 

§9. 

We shall now consider the more general transformation of nablas, viz. 
when the arguments are expressed in variable scalars and vectors which 
have certain geometrical relations between themselves. 

Take a to be a vector whose form is given by 

a = IX + m/j. + nv , (62) 

and write 

* = 's + "3S + 4' (68) 

where 

X/i = v , fxv = X , vX = p. . (64) 

Thus, let X, p, u be a rectangular system of unit vectors in which I, m, n are 
all functions of a position vector p which may be written in the form, 

p = ix + jy + kz , 

where i, j, k are as usual a rectangular unit system, being constant through all 
operations, although they are arbitrary at first choice. Since then p is the 
independent variable vector we may write at any time, 

9 . 9 . , 9 
$i +J 9l, + IC 9~z 



pV = * 5= +^ 3r . + *: 



= i i. + ,• A. + k £ 

dx J dy dz' 



(65) 



(67) 



> (69) 
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Then, if we take a linear vector function of the form 

nr = ^ISXt -\- ymS/jiT + \mSfT , (66) 

where v is used for p y when it is already applied ; we have the results 

pV = — ff <rV» ) 

„V = — t _1 pV- ) 
Laplace's operator becomes 

S p Vk„V = ^V^ftV . (68) 

where mt' = ttV is a self-conjugate linear vector function. 
To verify the results (67) we obtain, by using them, 

V 1 = — n*V l = — t (tVO = — ^ » I 

V m = — T<rV m = — T (o-V m ) = — ^ . 
VW = — W^l = — 7T („V W ) = — ^ > 

which are obtained at once from the definition of jr. 

Again, since we have the inverse function in the form 

— n~h = XS^nv^nr -f- /nSynyh + v£yZymr : Sylymyn , (70) 

we obtain by assuming the result of (67), 

„-V^ = — i-'pV^ = — *~\vl) = — * » 1 

^V^ = — T _1 P V m = — 7r -1 (V m ) = — /i , J> (71) 

„yw = — 7r-' P v» = — n-~ l (V n ) = — v . J 

which are also known from the definition of „.y. 

This transformation is useful in considering the orthogonal system of 
surfaces. Thus, if we have for these surfaces the equations 

F l { P ) = c l , F i {p) = c. i , F 3 ( P ) = c 3 , (73) 

then the conditions of orthogonality are 

S V F lV F 2 = , 
S V F N F S = , 
S V F iV F, = . 



> (74) 
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Since we may write 

V F 2 = m l x,\ (75) 



J 



the form of it is given by 

nz = V log TyFfivFs + V log T^F^F.r + V log TvF s S V F s r . (76) 

10. 

Perhaps it might be well to remark on the differential and variation sym- 
bols as they are usually expressed by quaternion notations. Thus, we write 

d = - Sda„ v , ) 

I (77) 

d = — Sdo„v . ) 

Here da and da are purely arbitrary vectors, they may take any values either 
in direction or in length, any variations in the values of da and da giving rise 
to corresponding variations in the values of d and d on the left side. 

For the sake of definiteness I propose to call d, as it is written in (77), an 
unconditional or total variation operator, and da an unconditional or total 
variation, and d and da, respectively, an unconditional or total differential 
operator and differential ; on the other hand, in 

d|« = — S8 u ojj , ) 

> (78)* 

<*i« = — Sd la a a y,) 

d la a and d la a are used to indicate that the vector a varies in a given direction 
parallel to a, and here the words conditional or partial are substituted for the 
words unconditional or total of the former cases. 

Since the conditional variation is the component of the unconditional 
variation in the direction of a, we have 

d u a = ar'Sada = aSa~ l da , ) 

> (79) 

d u a = a~ l Sada = aSa~ l da . ) 

Hence we obtain conditional variation and differential operators in terms of 
unconditional ones, 

d u = — Sar l a ^Sada = SVdaa' 1 Va„^ + d , ) 

> (80) 

d la — — Sa~\\?Sada = S Vdaa' 1 Fa„v + d . ) 

* Here I add the sign of parallelism, in order to distinguish the similar operator da in Tait's 
Quaternions, 3rd ed., p. 372. 
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Thus, in the case of orthogonal surfaces mentioned in the last section, writing 

vF t = mp = M, y (81) 

v F t = nv = N , J 
A, M, N are normals to the surfaces, and on account of orthogonality 

A\\d„p, M\\d ixP , N\\d l2 p, (82) 

where, for instance, d l2 p denotes the tangent to the curve of intersection of the 
surfaces F x and F 2 at the point p, or it is the change of the position vector in 
the direction of this curve of intersection. Hence we have to distinguish the 
three cases of the variation for each of the three normals according as the 
position vector changes along the curve 12, or 23, or 31 ; hence the importance 
of notation like (78). 
Thus, instead of 

dA — — A 1 Sdp p ^ l , 
we have 

d m A = dysA = — A x Sdntftfi = ~ 4&WpVi , 

= — A, SN -> pVl SNdp = A x 8 VdpN- 1 VN pVl + dA , 

d m A = d lA .A — — A^dtxptf! — — A l Sd iA p fl ^ l , 

H83j 
== — A.SA-'^SAdp = A^SVdpA-'VA^x + dA 

d m A = d u A = — A x Sd m p^ x = — A l 8d m p l ^ l , 

= — A x SM -» pVl SMdp = A X S VdpM~ l VM pVl +dA; J 

and similarly for M and N. 

Thus, for any function ii, we have 

dQ = d tA ,S -f d lVL £ + d lt ,S 

= d u Q + d u Q + d m SJ ; 
or, symbolically, 

d = d lk + <?,„ + d lr , (84) 

where X, p, v constitute any rectangular system of either constant or variable 
unit vectors, or not ; for 

d\K + d tll _ + d\* = — ^pV^ _1 Skdp — Sftfpr^Sfidp — S p yv~ l Svdp 

= — S^dp = d . 
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Thus the sum of conditional differentials or variations in three rectangular 
directions is equal to the unconditional differential or variation. 

Again, if we denote by pVn„ the component of p y in the direction of a, 
in the same manner that d^p has been used to denote the component of dp in 
the same direction, that is, if 

pVw = a ~ l &W = *<fo -I pV > ( 85 ) 

we may extend the definitions of (78), viz., 

#i<r = — Sdfrp P v = — Sdp p Vl<r . 



= — i>dp pVlcr , ) 

= — Sdp„^ w , S 



.. ( 86 ) 

die = — Sdfrptf 

where a may be any vector either constant or variable ; and here we must have 
the names unconditional and conditional nablas. 

But when the independent vector variable is required to move in a given 
direction t, we do not have necessarily 

d h = — Sdyptf = — Sd lr a^ . 

If we define Q to be a linear function of the argument C, 

Q{0 = Q(i) + QU) + Q(Jc) 

(87) 
where /, p., v form any variable unit rectangular system, then for any vector a, 



with 

Again, if 



dp = da 1 

^ = dp } (88) 



d ti = d . 



then 



da — tpdp , 
dp = f~ l da , 

trVllS = »V > pVlli = pV > 1 

<rV"f = ?'~\>vk . pVis — fWis ; }> ( 89 ) 

All these are old things in new dress, and so need no demonstration. 
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An unconditional variation, differential, or nabla is the maximum varia- 
tion, differential, or nabla, inasmuch as the conditional ones are components. 
Thus, if Q be the force, W the work done by the force, dp the element of the 
path, a the function of p being a given curve, then for the maximum work done, 
we have 

d W = - Siidp = - W.S^dp . 

However, when the path is constrained to be on a curve whose tangent is <r, 

d l<r W = - SQd uP = - Sli w d P , 

where Q w is the component of the force along the tangent of the curve a. If 
it happens that the force has a potential P, such that 

then 

dwW^-W.S^d^p, 

= — Siid la p = — PiS^d^p . 

Hence, whatever may be the curve a, we have 

W 2 -W 1 =P 2 -P l . 

As another illustration, we may prove the following theorem : — Any sys- 
tem of three surfaces which cut each other in their lines of curvatures, cut each 
other in constant angles. Let v„ v 2 , v 3 be normals to the three surfaces ; let 
d l2 p, d^p, d 3l p, be the tangents to the curves of intersection of the surfaces 1, 
2 ; 2, 3 ; 3, 1, respectively ; then the conditions that the surfaces 1, 2 cut each 
other in their lines of curvatures are 

= Sv^dv^d^ = S Vv x dv x Vv v v 2 = Sv l v 2 8v l dv 1 — v?Sv 2 dv x , 

= Sv 2 dv 2 d n% o = SVv 2 dv 2 FvjVj = Sv l v 2 Sv 2 dv 2 — v 2 8v x dv 2 ; 

whence we obtain 

Sv^dvi = — 2 Sv x dv x Sv x v 2 , 



Sv 1 dv 2 = — 2 Sv 2 dv 2 Sv ] v 2 , 



Adding these, we obtain 



Sv x v 2 ~ * 



dv? dv? 



or, 

d log S Vl u t = id log (v,V) = d log (T Vi n 2 ) 
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Integrating, we have 
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whence 

and similarly, 



Sv x v 2 = — c 3 Tv t v 2 . 
cos (v„ v 2 ) = c 3 , 
cos (v 2 , v 3 ) = C! , 
cos (v„ v t ) = c 2 . 



Of course this system of three surfaces includes the orthogonal system of 
surfaces. 

§11. 

For any function Q, linear in each of the arguments, McAulay* denned 
his very ingenious £-pairs in the form 

Q (C. = Q (,vi. ft) = Q (ft, pVO , 

and thus far I have used the pair in the same sense, but necessitated by the 
extension into the case of quaternions in general, I propose now to define 
£-pairs, jy-pairs, and £-pairs in the forms, 



Q (C. = Q UpVi. ft) = Q (£p» pVi) = etc., ] 

Q (?, r i)= Q (jaVi. ?i) = Q (Kq u Vi) = etc., 



(90) 



J 



(91) 



where ^ and ;o are any quaternion and vector, and McAulay's ("-pair is the 
negative of the same pair here. Employing a unit rectangular system, we have 

Q(C,0= -Q (*, »") - C (i, ./) - (*. *) . 

Q 0?, v) = Q (i. i) - Q (*. *") - <2 (i. i) - « (*, *) . 
£(£, €) = 0(i;i) + Q(i, i) + QUJ) + Q(k, k) . 

For the relation of these three pairs, we have 

Q(v> V) = 0(1, 1) + QC, = C(f, *5) , 

«(f, f) = «(i, i) - ec, o = c(7, *?) . 

«(7.7) + «(f,f) = 2G(l,l) . 

£0?, 9) -<?(?,?) = 2(? (CO - J 

* Utility of Quaternions, etc. , Arts. 1-4. 



(92) 
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Thus we obtain at once 

Si;r= 3, FCC = 0,^1 

Sw= 4 , V m =Q,\ (93) 

SH = -2, 7ff=.0.J 

The operators QS£(), rjSrj () operated on any vector and quaternion respec- 
tively produce the same vector and quaternion, while £S$ ( ) produces the 
conjugate ; thus, 

tjS W = p, I (94) 

$Ssp = Kp. , 

The operators £ Ff ( ), J? Fi? ( ), £ Ff ( ) produce twice the vector, thrice the 
vector part, and the conjugate of the vector part respectively ; thus, 

CFe<r = 2<r, 

7 Vrp = 3 F> , I (95) 

£F£p=^F> = - T>. 

The operators St^pSt^ ( ), SrjqSrj ( ), S$qS£ ( ) produce scalars of the product ; 

thus, 

S&Sro = Spa , 1 



SirjqSrjp = Sqp , 
SZqSZp = S^ . 



(96) 



The operators #£/> F£ ( ), etc., produce the scalars of the products of vector 
parts ; thus, 

8CpVZo = — 2Spo, 1 

SrjqVyjp = - SVqVp , j> (97) 

S?qV$p = SSVqVp.J 
The operators V£pS£ ( ), etc., produce vectors of the products ; thus, 
VZpSZa = Vap, I 

VrjqSrjp = — VKqKp =Vpq , [ (98) 

VSqSSp = - YKqp = VKpq . 
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The operator V£p VQ ( ) produces zero ; thus, 

VQp V^o = . 
The operators £f>S£ ( ) produce products ; thus, 

£pS£o —op, 1 

W8&=M, \ (99) 

ZgSsp — Kpq . J 

The operators 5^()^and VQ()Q, etc., produce upon any quaternion q and 
p, the following results : — 



(100) 



scq: = 


BSq, 


| 


Srjqtj = 


iSg, 


&im = o , > 

8$p6 = , J 


Ssq$ = - 


-2Sq; 


VQqr = - 


-Vq, 


PfrC = - /» , 


Vrjqrj = 


o, 


Vypy = , 


V$q$ = 


2Vq; 


Ffy? = 2p. 



(101) 



Cases in which more than one pair of C> )?, and $ occur are rather inter- 
esting when reduced to simpler cases ; thus, if Q is any function linear in each 
of the arguments, 

Q (C> Q ^GG = Q (Ci, Q Sc 2 C t = Q (C, , ] 



<?(ft> %) <%% = @0?i> %) %7i = Q(v> v) > 

«(fi, f.) «if, = £ (£., fj #£ = <2 (7, 7) • J 
If $ is a linear function in £,£ 2 , j^j, £j,£ 2 in combination, then 

0CiC.)FCi:, = — 2«(c) c, 

G(7i%)^i7a = °. 

e(Ci^7"«.= 2^(0 c, 

G(Wfc)7Wi = 40(C)C, 

0(^ 2 )Fc 2 ?i= 0. 



(102) 



(103) 
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Thus, if Q be linear in the combination CtfCz; where q is any quaternion, 

Q (dsG ^:,r, = (:*&) <%g = 30 m - ( r ? ) , i 

(Ml*) &W* = (M%) %=fi = 4 (^) . f ( 104 ) 

(%£ 2 ) -SfA = £ 62Q 5f A - 40 (ty) . J 

«(Ci?C») FCCi = 20(1) Fy + 20(OCty , 1 

(M%) ^?i% = 40 (1) Fy + 40 (C) Cty • ^ (105) 

I 

<?(£#£,) F?A = o J 

(CtfG) FC2C. = - 20 (1) Vq - 20 (0 t$q , 1 

«(M%)^i% = 0. > (106) 

(%£ 2 ) F?A = - iQ (1) F? - 40 (0 C<% • J 

The definitions of two pairs of £, j?, £ can be at once inferred from those 
of one pair, and worked thus far on that assumption, but they are written in 
the forms for any function Q, linear in each of the arguments, 

0(Cl. <Tl» C* Q = 0(oVl » Pl> <rV2 . ^2) = Q(P\, pVl » ^2, »V*) = 6tc -» 1 

I 

(Vl Vu %, %) = (i-«Vi, ?i, 1&V2, ^2) = (?i, jraVi. A. KpVi) = etc-, ^ (107) 

0(?i» fi» f», ?*) = 0(«Vi , ?., pV2. />*) = 0(?i, ,Vi >i>2> PV2) = etc., J 

where /) and a are any vectors, which may be the same or different ; p and q 
may also be any quaternions. 

Any self-conjugate linear functions may be written 

fP = "{Sep , ) 

> (108) 

0q = a^S-rtf , ) 

where the suffix of the scalar a is supposed to vary with Q and rj. 
Any linear function may be written in the form 



<PP = f'SCP . ) 
0q = (ptjSrjq 



V (109)* 



* Hamilton's Element's, Art. 365. 



KIMURA. ON THE NABLA OF QUATERNIONS. 

Their conjugate functions may be written 

f'P = Z S P<PZ > ' 
(p'q = rjSqtir) 
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mo) 



The linear vector function and its conjugate derived from a linear quaternion 
function may be written 



9P = VVQSrp , \ 
tp'p = rSp<Pr . S 



(111) 



The invariants and rotation vectors of <p (and of its Hamiltonian inverse func- 
tion) may be written, 



Thus, if 



m' = 8ZK , 

m" = S£<p£ ; 

If <p and it are any two linear vector functions, 

Stp&tQ = Siztftfpi = StptfMi = SpVifapi = Stfi^fpi ' 
= <SjOif'TpVi = Sp\rf<Ptfi 

v<pzk = yipwpi — — vxtfrfPi 

= VfPiXpVi = — VnpifpVi- 

Q (?» 2> ?> 9> • • •) 



(112) 



(113) 



be any quaternion function linear and homogeneous in q, of the wth degree, 
then 



— , (SpxtViY Qi is, q,q,q,---) = Q{p, p> p,p,---), 



(114) 



where p is any quaternion constant with respect to q, since 

(8px«Vi) Qi (a, q,q,---) = Q (ySpy, q,q,---) + Q (a, ySpy, q,---) 

+ Q (q, q, ySpy, ?>•••)+ • • • (» terms). 
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When these are operated on again by Sp KlJ y each term produces n — 1 terms, 
and so on, until finally after the nth operation, all ^'s within Q are changed into 
^>'s ; hence the result. Thus, for the function Q, which is linear in q or p, or 
homogeneous of the rath degree in q or />, we obtain, if a and a are respectively 
a constant vector and a constant quaternion, 



1 



(-l) n ^(Sa^rQ<P) =Q(«), 



i-ASa KqV TQ{q)^Q{a) ) > 



±,(S**rQ(2) =Q{Ka). 



(115) 



This is a generalization of a problem in Tait's Quaternions.* 

Hamiltont had proved that if dq is a constant quaternion different from 
zero such that d i q = 0, d s q = 0, . . . , then 

e"f(g)=f(g + dq), 

where y is any function, provided the first m differentials do not become infinite. 
But since the operator Sp Kq ^, in which p is a constant with respect to q, is 
equivalent to the operator d with the condition added that dq = the constant 
p, we may write in our notation the Taylor's Theorems in the forms : 



e s P^f(q)=f{q+p), 1 
e^f{q)=f{q + Kp), Y 



(116) 



Similarly 

(<**»** = €**«) f(q)F(j>) =f( q + KpVi )F{p,) = f(q 1 )F(p + KtVx ), 1 

(«*w = e s«™)f(q)F(p) =f(q + pVl ) F{ Pi ) = /(<?,) F(p + , Vl ), [ (117) 

^V/(/>) F{o)* =/{p + , Vl ) F(a,) =f(p,) F{a + pV i), J 

where p and q, a and p are independent of each other ; and by changing the 
signs of p, <r, xjtf, p y, <rV i Q b°th members we have the differences instead of 
sums in the arguments of the functions f and F. 

*3rded., pp. 420, 400, 399. 
+ Elements, Art. 342. 
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In the last section we gave the notation 

di* = 8d la p upV = Sdp KpVu , 
— 8a jrpV 8a~ x dp = . . . , 

to mean that it gives by operation on a function of p the differential of the 
same function when the variable p moves parallel to a given vector a ; if the 
vector dp is not only parallel to a, but equal in tensor to the same, then 
the above becomes Sa Kl ^, whatever length the vector a may have. But, 
Sa-KpV i g 8dp Kf ^ = d, in which the arbitrary vector dp is changed to a con- 
stant vector a ; hence the perfect accordance of both notations. Here it might 
be remarked also, that if 8a K rf is operated only once, then a may also be 
any function of p. 

New Haven, October 25, 1895. 



